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ABSTRACT 

The effect of massive neutrinos on the evolution of the early type galaxies 
in size and stellar mass is explored by tracing the merging history of galaxy 
progenitors with the help of the robust semi-analytic prescriptions. We show 
that as the presence of massive neutrinos plays a role of enhancing the mean 
merger rate per halo as well as the merger driven increment in halo mass, the 
high- z progenitors of a massive descendant galactic halo evolve more rapidly in 
mass-normalized size for a AMDM (A Cold Dark Matter + massive neutrinos) 
model than for the ACDM (A Cold Dark Matter) case. We provide a physical 
reason for why the halo mass growth rate and the merger rate are higher in a 
AMDM cosmology and conclude that if the presence and role of massive neutrinos 
are properly taken into account it may explain away the anomalous compactness 
of the high- z massive ETGs compared with the local giant ellipticals with similar 
stellar masses. 

Subject headings: cosmologyitheory — large-scale structure of universe 



1. INTRODUCTION 

It has been recently discovered that the quiescent early type galaxies (ETGs) at high red- 
shifts (z > 1) are smaller in size (i L) by a factor of (3-5) compared with the local giant ellipti- 



cals o f similar stellar masses (M*) (IDaddi et al.ll2005l ; iTrujillo et al.ll2006l ; Ivan Dokkum et al. 



20081 ). which implies that the high- z massive ETGs must have evolved very strongly in size 
but not in stellar masses. In the beginning it was suspected that the stellar masses of the 
high-;? ETGs from the photometric data should be overestimated. But the subsequent spec- 
troscopic estimates of their dynamical masses disproved this suspicion, confirming the rapid 
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size growth of the high-z massive ETGs (Ivan de Sande et al.ll201ll ). It has been a puzzling 
mystery what mechanism drove such strong size evolution of the ETGs without triggering 
their star formation activity. 



Fan et al.l (120081 120101 ) proposed a scenario that the compact sizes of the high- z massive 
ETGs should be attributed to the AGN feedback effect of blowing off the baryonic gas out 
of the halo potential wells. The merit of this scenario was that it can explain naturally 
why the compact E TGs are observed at high redshifts z > 1 when the quasars were most 
energetic. However, iRagone-Figueroa fc Granatol (1201 ll ) cast a doubt on this scenario based 
on the argument that it is difficult to understand in this scenario why the most compact 
high-z ETGs are the quiescent ones possessing old stellar populations. 

An alternative popular scenario was that the progenitors of the present local giant 
ellipticals underwent very fr equent dry mergers through which they grew rapidly in size but 



only m il dly in stellar mass (IHopkins et al. 



2009 



Khochfar & Silk! 12009 : iNipoti et al.l 12003 



2009al lbl; lOser et aL 2012 ; van der Wei et al. 2009 ). At first glance this scenario seems quite 
plausible, fitting well into the standard theory of hierarchical s tructure formation based 
on the ACDM (A + cold dark matter) cosmology. But, recently INipoti et al.l (120121 ) have 
revealed that this scenario in fact challenges the standard model rather than supporting it 
since even the maximum merger-driven size growth predicted by the ACDM cosm ology is not 
fast enough to catch up with the observational trend (see also lCimatti et al.ll2012f). unless the 



key c osmological parameters deviate substantially from the WMAP7 values (IKomatsu et al. 
201 lh . 



We note, however, that when INipoti et al.l ( 120121 ) calculated the maximum merger- 
driven size growth consistent with ACDM cosmology merger rate, they missed one cru- 
cial point that the real structure formation in the Universe may not proceed in a per- 
fectly hierarchical way due to the presence of small amount of hot dark matter compo- 
nents — massive neutrinos. Since it was found and confirmed that the neutrinos are not 



massless (Ahmad et al. 



2001 



Arnett fc Rosnedl 19871 : ICleveland et al.lll998l : iFogli et al.ll2006 



Fukuda et al.lll998l : iMaltoni et al.l 120(5 )7 plenty of theoretical works have been devoted to 



investigating what effect the massive neutrinos would have on the formation and evolution 
of the cosmic structures an d how significant the effect would be (for a recent review, see 
Lesgourgues fc Pastorll2006l ). 



It is now understood that on the time-dependent scales where the massive neutrinos 
become non-relativistic, the high velocity dispersion of massive neutrinos would damp out 



DM halos ( 


Abazaiian et al. 


2005; 


Agarwal & Fe 


dman 


2011: 


ArhiDova et al. 


2002; 


Bird 


et al. 


2012; 


Ichiki & Takada 


2012 


Lesexmreues et al. 


2009 




Mantz et al. 


2010; 


Marulli et al. 


2011; 
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Saito et al.l I2008L 120091 ). The recent works of ISong fc Led (1201 ll . 120121 ) also suggested that 
the effect of massive neutrinos should also bring about non-negligible change in the mean 
merging rate of DM halos. Assuming a strong dependence of the size growth of the observed 
ETGs on the merging rate of their underlying DM halos, it is logically expected that the 
ETG compactness would be also affected by the presence of massive neutrinos. 

Here we claim that the presence of massive neutrinos (z/) could significantly promote 
the size growth of the massive ETGs, explaining away the observed anomalous compactness 
of massive ETGs at high redshifts. To prove our claim , we trace backwar d the progenitor 
evolution of a massive descendant galactic halo, just as iNipoti et al.l (120121 ) did, but taking 
the effect of massive neutrinos into account. The organization of t his Paper is as foll ows. 
In section [2] we provide a brief review of the model constructed by INipoti et al.l ( 120121 ) for 
the size and stellar mass evolution of the ETGs. In section [3] we present a new model for 
the size and stellar mass evolution of the ETGs for a AMDM (A Mixed Dark Matter which 
represents CDM + v) cosmology and explain how the observed compactness of the massive 
high- z ETGs can be explained away by taking into account the presence of massive neutrinos. 
In section @] we summarize our result and discuss its caveat and cosmological implication as 
well. 



2. MODEL OF THE ETG PROGENITOR EVOLUTION: A REVIEW 



Nipoti et al.l ( 120121 . hereafter N12) envisaged a simple picture where a dissipationless 
dry merger (which maximizes the size growth of the progenitor galaxies) occurs between a 
main galaxy of mass M main (z) and a satellite galaxy of mass M sat (z) to form a merged galaxy 
of mass Mh(z) at redshift z, which represents a progenitor of a descendant galaxy of mass 
Mh(zo) observed at a given redshift zq < z. Then, they modeled the merger-driven evolution 
of the size and stellar mass of the progenitor galaxy as 



d 2 M+ dM, 



(d 2 M h \ 



merg 



dzd£ dM h \ dzd^ J ] 
d 2 \nR e d\nR e 1 dM± fd 2 M h \ 
dzdi ~ d In Af* M* dM h \ dzd£ J 



(1) 



(2) 



merg 



where £ = M sat /M main and the subscript "merg" in the right-hand side means the change 
caused by the halo merging events not by the smooth accretion of matter. Here the progenitor 
size R e represents the effective radius that encloses half the progenitor stellar mass. 

To evaluate Equations ([I])-©, N12 prescribed the three ingredients: (i) the evolution 
of halo mass d 2 M} l j[dzd^)\ (ii) the stellar-to-halo mass relation dM^/dM^, (iii) the stellar 
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mass to size relation dM*/dR e . The first ingredient d 2 Mh/(dzdC,) was prescribed as 

fd 2 M h ^ iM h {z) d 2 N n 



< merg 



\ dzdt; J 



merg 



1 + £ dzd£ 



(M h ,z,0 



(3) 



where £Mh(z)/(l + £,) is the amount of mass that a progenitor gathers during a merging 
event (i.e., the merger-driven increment in halo mass) that occurs in a redshift interval of 
[z, z + dz] with satellite-to-main mass ratio £, while d 2 N merg / (dzd^) represents the mean 
number of such merging events per a merged halo of mass Mh(z) (i.e., mean merger rate per 
halo). From here on, we denote the mean merger rate per halo by B(Mh, z,£). 

The tricky part in evaluating Equation ([3]) is the quantity Mh(z) in the right hand side, 
which represents the halo accretion history, i.e., the total halo mass of a progenitor at z, which 
grows not only through merging between the main and the satellite galaxies but also through 
diffusive accretion of matter. Thus, the functional form of Mh(z) has to be first determined 
for the evaluation of Equation fl3]) rather than regarding it as a integration va r iable. N12 used 
the fitting formula for B and Mh(z) derived numerically bv lFakhouri fc Ma (2008, hereafter 
FM08) from the Millennium simulations for a ACDM cosmology (jSpringel et al.ll2005l ). 



As for the prescription of the second ingredient dM^/dM^, admitting that no standard 
model has yet to be established for t he stellar-to-hal o mass relation (SHMR ), N12 considered 
three different SHMRs pro vided by iBehroozi et all ( bold . BIO hereafter), iLeauthaud et al.l 



(120121 . L12 hereafter) and IWake et al.l (1201 ll . Wll hereafter), respectively. These three 
SHMRs were all obtained by matching various physical and statistical properties of the 
galactic halos from the high-resolution simulations to those of the observed galaxies. The 
difference among the three SHMRs comes from the different redshift ranges to which the 
matching between the numerical and the observational results was applied. In the current 
work, we consider only the first two SHMRs (BIO and L12), excluding the last one (Wll) 
which is valid only for relatively narrow redshift range. 

The BIO and the L12 SHMRs have the following same functional form in the redshift 
range of 1 < z < 2: 



logil4(i\4) = logMi + glog 



M^ ) 1 + (M*/il4, )~ 



(4) 



Here the values of the five model parameters, Mi, M* , q, P, and 7, vary with the scale 
factor a as 



logMi(a) 
logM* j0 (a) 
q(a) 



Mi,o + M 1>a (a-l), 
M*, ,o + M* Aa (a - 1) 

qo + q a (a - 1), 



(5) 
(6) 
(7) 
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p(a) = Po+Pa(a-l), 
7(a) = 7o + 7a(a-l). 

At z > 1, M* )0 becomes time dependent as 

logM^o(a) = M* j0 ,o + M^ a (a - 1) + M* Aa2 (a - 0.5) 2 . 



(9) 



(10) 



The best-fit values of these model parameters are provided in Table 2 in BIO and Table 5 in 
L12, resp ectively. Although the original best-fit parameters of the L12 SHMR is valid only 



Nipoti et al.l ( 120121 ) extrapolated the validity of the L12 SHMR to higher redshifts 



at z < 1, 

of z > 1 and redetermined the best-fit parameters. 

The third ingredient, the size-to-stellar mass relation dR e /dM+, was prescribed as 



dhiR e 
din 



2 - 



in(i + e 1 - 4 ; 

ln(l + 



(Naab et al. 


2009; 


Oser et al. 


2012) 



shapes; the satellites follow the parabolic orbits; there is no energy loss during the dissipa- 
tionless dry merging events. Plugging the prescribed three ingredients into Equations ([I])-© 
and integrating them over £ and z, N12 have finally made the quantitative predictions of the 
standard ACDM cosmology for the merger driven evolution of the size and stellar mass of 
the progenitors R e (z) and M*(z) of a given descendant halo at zq < z. Plotting the locations 
of the high-z progenitors at z > 2 in the R e (z)-M*(z) plane and comparing them with the 
local relations (see Figure 10 in N12), they demonstrated clearly that the sizes of the high-z 
progenitors in ACDM universe is much larger than the observed ones. 

Estimating the uncertainties in R e (z) and M+(z) due to the simplified assumptions that 
they made to prescribe dR e /dM+ and dMh/dM+, N12 showed that the most severe uncertain- 
ties come from the SHMRs and that the disagreement between the ACDM prediction and 
the observational result is robust against the SHMR uncertainties. In other words, it was 
confirmed that although several simplified assumptions were made in their determination of 
the size and stellar mass evolution of the progenitor galaxies, the large discrepancy between 
the ACDM prediction and the observational result is not due to the inaccurate modeling of 
dRe/dM* and dM h /dM* . 



In the follow-up work, ICimatti et al.l (120121 ) quantified the compactness of a progeni- 
tor galaxy as the mass-normalized size of R e M^ with M u = M*/(lO n M ), and mod- 
eled the evolution of the progenitor compactness as a power-law scaling with redshift: 
i? e M~°- 55 oc (1 + zf, given that the effective s izes of the observed g alaxies scale as apower- 



law of their stellar masses, R e oc M®' 55 (see also lNewman et al.ll2012l . and references therein). 
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Cimatti et al.l ( 120121 ) determined the best-fit slope to be ~ —0.6 for the ETGs in the red- 



shift range of < z < 2 and to be of ~ —1 for the ETGs at < z < 2.6, confirming the 
rapid evolution of the ETG compactness at high redshifts. 



EFFECT OF MASSIVE NEUTRINOS ON THE ETG PROGENITOR 

EVOLUTION 



In this section we investigate how the presence of massive neutrinos affects the com- 
pactness evolution of the massive ETGs, applying the N12 model to the case of a AMDM 
cosmology. Basically, we compute Equation (JT|)-([2j) in a similar manner but with a modified 
prescription for the first ingredient, d?Mh/dzd£, by incorporating the presence of massive 
neutrinos into the picture. First of all, we would like to determine the mean merger rate per 
halo B and the halo accretion history Mh(z) for the AMDM case. Of course, the most accu- 
rate way to determine these quantities should be to run repeatedly high-resolution N-body 
simulations for AMDM models with various different values of the neutrino mass ^2 m v and 
then to determine the empirical fitting formula for B(z) and Mh(z) as a function of 
by analyzing the N-body results. Since such simulations are not available at the moment, 
we utilize the less accurate but practical semi-analytic formula for B and Mh derived from 
previous works. 



Zhang et al.l (120081 ) analytically derived the following formula for the mean merger rate 



per halo by inc orporating the ellipsoida l collapse dynamics i nto the extended Press-Schechter 
(EPS) theory jPress fc Schechterill97i [Lacev fc Colelll993h : 



B(M h ,z,g) = B Bph (M h ,z,£) x A exp - 



1 + A 2 S* /2 



1 



r(3/2) 



(12) 



where A n = 0.866(1-0. 133u r 



-0.615\ 



Ai = 0.308za 



-0.115 



A, 



0.0373^ -°- 115 , v 



u 2 {z)/S{M h ) 



and Si = ASj/ S( M h ). Here, B Kp h repre sents the original EPS model based on the spherical 
collapse dynamics (ILacey fc Coldll993l ): 

dd c (z) Ml dS(Mi) 1 



B sph (M h ,z,0 



where S c (z) = 8 C /D(z), S c 



dz {l + ifMi dMi ASi^TiASi) 1 / 2 ' 
1.68, D(z) is the linear growth factor, ASi = S(Mi) 



(13) 
S(M h ), 



and S(M) = a 2 (M) is the variance of the linear density field smoothed on the mass scale of 
M. The quantity Mi in Equat ion (fT2l repr e sents the mass of a merging halo which could 
be either a main or a satellite. IZhang et al.l ( 120081 ) showed that Equation ( JT2|) significantly 
improves its spherical counterpart, Equation ( JT3|) . agreeing much better with the N-body 
results. 
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We adopt this ellipsoidal EPS model to evaluate the mean merging rate per halo B for 
a AMDM cosmology whose linear density power spectrum is characterized by the neutrino 
mass fraction, f„ = Q v /Q m (where fl m and fl u are the matter and the massive neutrino 
density parameter, respectively). Extrapolating the validity of the ellips oidal EPS model to 
AMDM cosmology may be justified by the results of iMarulli et al.l (120111 ) which showed that 
the PS-like approaches are valid for the evaluation of the halo mass function even for AMDM 
cosmologies. For the AMD M linear density power s pectrum and growth factor, we utilize the 
analytic formula given by lEisenstein &: Hul ( 119991 ) . To normalize the AMDM linear power 



spectrum, we make its amplitude satisfy as = 0.8 for the case of f v = (i.e., without massive 
neutrinos) , which is equivalent to the large -scale normalization. The other key cosmological 
parameters are set at the WMAP7 values ( IKomatsu et al.l 1201 if ). 



Figure [T] shows the mean merger rate per halo as a function of £, for three different 
cases of f v with M h (0) = 10 13 M o (which corresponds to M*(0) « 10 u /i -1 M o ). As can be 
seen, the larger the neutrino mass fraction is, the higher the mean merging rate per halo is. 
For the case of f v = 0.05 (which corresponds to the WMAP7 upper limit of ra v < 0.58eV 
with effective number of massive neutrino N e ff = 4.34), the overall merger rate per halo 
increases by a factor of 3 relative to the case of f v = 0. It is worth mentioning here that the 
redshift z at which the descendant galactic halo is observed is set at the present epoch, i.e., 
zq = 0, unlike zq = 1 in the original N12 model, as we are interested in the full progenitor 
history from the high redshifts to the present epoch rather than focusing on explaining the 
high-z phenomena. 

As for t he mass acc r etion history for a AMDM universe, we use the analytic formula 
proposed by IZhao et al.l (120091 ) . Noting the existence of a simple expression f or the mass 



accret ion rate as a function of halo mass, redshift and cosmological parameters, IZhao et al. 



(120091 ) developed a theoretical model for the mass accretion history which has the following 
universal form: 

d\oga{M) u{z, M) - p{z, z ohs , M h {z )) 



d\ogS c (z) 



5.85 



Here, 



u(z,M) 
p(z, z , M h ) 
p(z ,M h (z )) 



S c (z) 



o-(M)10 dl °s <J / dl °s M ' 
p(z Q ,M h (z Q )) x Max 
1 



0.272/a;(^o,M h (z )) 
u[z , M h (z }) 



(14) 

(15) 
(16) 
(17) 



l + [Lo(z ,M h (z ))/4] e 2 

Note that in Equation (j!4p . l/a(M) is used as a mass-like variable and l/5 c (z) is used as a 
time-like variable. 
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Since IZhao et al.l (120091 ) showed that Equations flT4l)- f|T7j) works on wide mass range 
at various redshifts for several different cosmological models, the applicability of Equations 
( TT1]) -(1T7 ]) to AMDM cosmology is expected. Expressing Equation (THj) in terms of mass and 
redshift according to the chain rule as, 

dM h dM h d\og5 c (z) dloga(M h ) 



dz dloga(Mh) dz d\og5 c (z) 



and integrating Equation (1181) over z, we finally determine the halo mass accretion history 
Mh(z) for a AMDM universe. Figure [2] shows the evolution of and M* via dry mergers 
with £ > £ min = 0.03 for three different cases of f u in the left and right panels, respectively. 
The mass of a descendant halo at present epoch is set at Mh(zo = 0) = 1O 13 M (corre- 
sponding to ~ 10 11 Mq) for all three cases. Here, we follow the evolution up to redshift 
z ~ 1.5 since the formation epoch of a descendant halo with mass M/ l (z = 0) = 1O 13 M 
is approximately z = 1.5 . As can be seen, both of Mh and M* evolve more rapidly for the 
case of higher value of f u , and the differences among the three cases of f v in Mh and M* 
become larger at higher redshifts. This result indicates that the high-z progenitor galaxies 
of a descendant halo with the same mass accrete larger amount of DM and stellar masses 
via dry mergers in a AMDM universe than in the ACDM case. 

The physical reason why the halo mass growth rate and the merger rate are higher 
in a AMDM cosm ology may be related to the later formation epochs of dark halos (e.g., 



Song fc Ledl201ll ). The dark halos tend to form later in a AMDM universe than in a ACDM 
universe since the small-scale powers are reduced due to the free streaming effect of massive 
neutrinos. In a ACDM universe the small galactic halos form much earlier than the larger 
halos and thus the interval between the formation epochs of small galactic halos and larger 
halos is long. However, in a AMDM universe the formations of dark halos (small and large 
halos alike) are delayed and thus the interval between the formation epochs of small and 
large halos is much shorter than in the ACDM case. In other words, the small halos need to 
merge faster into the larger halos in a AMDM universe. 

Now that we have determined the mean merger rate per halo and the mass accretion 
history for the AMDM case, we plug them into Equation ([3]) to prescribe the first ingredient 
of our model. Regarding the second and third ingredients (i.e., dR e /dM+ and dM^/dMh), 
we use the same prescriptions i.e., Equations (jll)- (fTTT) . of the original N12 model under the 
assumption that these relations are still valid for the AMDM case. It is worth, however, 
discussing whether or not this assumption is reasonable here. As mentioned in section [2j 
the SHMRs were all determined by matching the observational data to the numerical results 
by applying several different techniques such as abundance matching, gravitational lensing 
and so on, some of which are not based on ACDM models. Therefore, even though the 
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SHMRs that N12 adopted were obtained for the ACDM case, we expect it to work for 
the AMDM case. Regarding the size-to-stellar mass relation, Equation ( II ip . it is a local 
relation obtained by applying the first order galactic dynamics which is independent of the 
background cosmology as long as the gravity is well described by the Newtonian dynamics. 
Therefore, it should be also valid for the AMDM case. 

Plugging the prescribed ingredients into Equations (Jl])-(j2j), we finally calculate the size 
and stellar mass evolution of the progenitor galaxies of a descendant halo observed at present 
epoch in a AMDM universe and then determine the compactness of the progenitor galaxies 
at each redshift as R e M^' 55 . Figure [3] shows how the compactness of the progenitor galaxies 
of a descendant halo observed at z = evolves in a AMDM universe. The left and right panel 
corresponds to the cases of M h (0) = 10 13 h~ l M Q and M h (0) = 5 x 10 12 h~ l M e , respectively. 
In each panel, the solid, dashed and dot-dashed lines correspond to the cases of f u = 0, 0.02 
and 0.05, respectively. As can be seen, the high-2; progenitors become more compact for 
the higher-/,, case, which confirms that the presence of massive neutrinos plays a role of 
promoting the size growth of the massive ETGs at high redshifts. Furthermore, this effect 
of massive neutrinos on the size evolution is stronger for the case that M/ l (0) has a higher 



massive ones ( 


Cimatti et al. 


2012; 


Daddi et al. 


2005; 


NiDoti et al. 


2012 


van Dokkum et al. 


2008 





The results shown in Figure [3] have been obtained by using the B10 SHMR. To examine 
whether or not the use of a different SHMR will change the trend, we repeat the calculation 
of the compactness evolution but with using the L12 SHMR. The top panel of Figure HJshows 
the compactness evolution obtained by using the L12 SHMR as thick solid and dashed lines 
for the cases of f u = and /„ = 0.05, respectively. The results shown in the left panel 
of Figure [3] (only for the two cases of f v = and 0.05) are also plotted as thin lines for 
comparison. Since the L12 SHMR was originally obtained for the galaxies at z < 1, we 
show the results only at z < 1. The bottom panel of Figure H] shows the fractional difference 
between the two SHMR cases. As one can see, the use of a different SHMR does not destroy 
the trend that the progenitor compactness evolves faster for a higher f v case, which indicates 
that our result is qualitatively robust against the uncertainties of SHMR. However, it has 
to be also noted that using a different SHMR yields quantitatively a different compactness 
evolution of the progenitors and thus it will be important to refine the SHMR as accurately 
as possible. 

We have so far studied theoretically the progeni tor evolution of one single descendant 
halo existent at z — . The observed ETGs shown in lCimatti et al.l (120121 ) must correspond 
not to the progenitors of one single descendant halo existent at z = but to the progenitors 
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of different descendants existent at different redshifts. Therefore, to make a comparison with 
the observational results, it is necessary to model theoretically the progenitor evolutions of 
different descendant halos existent at different redshifts for each case of f v . We consider only 
those progenitors w hose distributi o n in t he M±-z plane are similar to the observed ETGs 
shown in Figure of ICimatti et al.l (120121 ). having stellar masses in the range of 10 10 ' 5 < 
M*/M Q < 10 12 at redshifts of < z < 2.5. It amounts to considering those halos as 
representative descendants whose total mass lie between 5 x 10 12 M Q and 2 x 10 13 M . 

Approximating the compactness evolution of the theoretically modeled progenitors as 



the scaling relation of R e M x l 



0.55 



'1 + z) 13 , we determine the best-fit value of (3 with 



the help of the \ 2 minimization method for each case of f u . Figure [5] shows the best-fit 
scaling relations (solid line) and compare them with the average compactness evolution of 
the modeled progenitor galaxies (dots) for three different cases of f u . As can be seen, the 
absolute value of (3 i ncreases from 0.89 to 1.02 as f v increases from to 0.05. According to 
Cimatti et al.l (120121 ). the absolute value of (3 obtained from the observed ETGs at redshifts 
of < z < 2.5 is close to unity. At face value it suggests that the observational result is 
consistent with the theoretically estimated value of (3 for the AMDM model with f v = 0.05. 



4. DISCUSSION AND CONCLUSION 

We have shown that the presence of massive neutrinos plays a role of enhancing the 
mean merging rate per halo as well as the merger-driven increment in halo mass on the 
massive galaxy scale and thus that the ETG compactness evolves much more rapidly in a 
AMDM universe. It is worth discussing here whether the AMDM prediction of the higher 
merging rate per halo is consistent with the observational indication. N12 mentioned that 
since the inferred merging rate from the observations of binary galaxy systems is even lower 
than the ACDM prediction, enhancing the merging rate should not be the key to explaining 
away the anomalous compactness of the high-z ETGs. However, very recently, iJian et al. 



( 120121 ) have demonstrated that the merging rate inferred from the observed galaxy pairs 



might be severely contaminated by the spurious projection effect. 

However, our current model is only a reasonable approximation based on several sim- 
plified assumptions to the true effect of massive neutrinos on the ETG evolution. It will 
be definitely necessary to improve and refine our model by incorporating more realistic as- 
sumptions and adopting more accurate prescriptions. In the first place, it will be essential to 
determine more accurately the mean merger rate per halo by using high-resolution N-body 
simulations for a AMDM universe. Although the ellipsoidal EPS model that we have used 
to calculate the mean merger rate per halo is a significantly improved version of the original 



spherical EPS model, it still suff ers from maximum 20% errors when compared with the 



N-body results ( jZhang et al.ll2008l ). The size-to-stellar mass relation also has to be improved 
by accounting for the following realistic aspects of the true merging process of progenitor 
galaxies: (i) the occurrence of dissipational wet mergers; (ii) the time lapse between the 
moment of halo merging and that of galaxy merging; and (hi) the presence of disc shaped 
progenitors. High-resolution gas simulations with massive neutrinos included will be required 
to address these complicated issues. 

Another thing that may deserve discussing is the question of whether or not the AMDM 
model is the only solution to the anomalous strong size growth of the high-z massive ETGs. 
One might think that different cosmologies such as models with primordial non-Gaussianity, 
dynamic dark energy or modified gravity might also affect the mean merger rate per halo, thus 
influencing size growth of the massive high-z ETGs. To examine if any other cosmological 
models could be also a solution, however, it will require much more works, totally renewing 
all prescriptions since in these models the extended EPS formalism is no longer valid, the 
gravitational dynamics changes and the SHMR could be also quite different. 

Our final conclusion is that the tension of the standard structure formation scenario with 
the observed anomalous compactness of the massive high-z ETGs can be greatly alleviated 
without changing the key cosmological parameters from the WMAP7 values if the presence 
and effect of massive neutrinos are properly taken into account. 
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Fig. 1. — Mean merger rate per halo S versus the mass ratio £ for three different cases of 
the neutrino mass fraction, f u . 



- 16 - 




Fig. 2. — Evolution of the halo mass (left panel) and the stellar mass (right panel) accreted 
from redshift z only through mergers with mass ratio £ > £ m i n = 0.03 for three different 
cases of the neutrino mass fraction: (/„ = 0, 0.02 and 0.05 as solid, dashed and dot-dashed 
lines, respectively). 
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Fig. 3. — Compactness evolution of the progenitor galaxies for the case of the massive 
descendant halo with M^(0) = 10 13 M Q (left panel) and for the case of the less massive 
descendant halo with M h (0) = 5 x 10 12 M Q (right panel). In each panel, the solid, dashed 
and dot-dashed lines correspond to /„ = 0, 0.02 and 0.05, respectively. 
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Fig. 4. — (Top panel): Comparison between the compactness evolution results obtained by 
using two different SHMRs: the case of using BIO SHMR as thin lines while the case of using 
L12 SHMR as thick lines. For each case, the solid and dashed lines correspond to /„ = and 
/„ = 0.05, respectively. (Bottom panel): Fractional difference between the thin and thick 
curves shown in the top panel for the two cases of f u . 



- 19 - 



10 



u 

Q_ 

ID 

d «- 



1 ■ 



f„=0.00 


f„=0.02 


f„=0.05 : 




^^^^ 


* •v^.v.* • 

— -VT 


• • 

: /? = -0.89 

.... i .... i .... i .... i .... i 


: /? = -0.93 

.... i .... i .... i .... i .... i 


: /? = - 1 .02 "v 



0.00.51.01.52.02.5 0.5 1.0 1.5 2.0 2.5 0.5 1.0 1.5 2.0 2.5 
z z z 



Fig. 5. — Average compactness of the progenitor galaxies (dots) with stellar masses in the 
range of 10 105 < M+/M & < 10 12 for three different cases of f u . The average is taken over 
the progenitor histories of representative descendant halos whose progenitors have the same 
stellar mass distribut ions in the redshift range of < z < 2.5 as the observational results 
shown in Figure 2 of ICimatti et al.l (120121 ) . In each panel, the solid line corresponds to the 



best-fit scaling relation of Re/M^ 
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